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INVARIANCE PRINCIPLES FOR OPERATOR-SCALING GAUSSIAN 

RANDOM FIELDS 

HERMINE BIERME, OLIVIER DURIEU, AND YIZAO WANG 


Abstract. Recently, Hammond and Sheffield fl7] introduced a model of correlated random 
walks that scale to fractional Brownian motions with long-range dependence. In this paper, 
we consider a natural generalization of this model to dimension d > 2. We define a Z d - 
indexed random field with dependence relations governed by an underlying random graph 
with vertices Z d , and we study the scaling limits of the partial sums of the random field over 
rectangular sets. An interesting phenomenon appears: depending on how fast the rectangular 
sets increase along different directions, different random fields arise in the limit. In particular, 
there is a critical regime where the limit random field is operator-scaling and inherits the 
full dependence structure of the discrete model, whereas in other regimes the limit random 
fields have at least one direction that has either invariant or independent increments, no 
longer reflecting the dependence structure in the discrete model. The limit random fields 
form a general class of operator-scaling Gaussian random fields. Their increments and path 
properties are investigated. 


1. Introduction 


Self-similar processes are important in probability theory because of their connections with 
limit theorems and their intensive use in modeling, see for example [44] . These are processes 
(X t )te r that satisfy, for some H > 0, 

(1) (X(Xt)) t£R f = A H (X(t)) tm , for all A > 0, 

where ‘ = ' stands for ‘equal in finite-dimensional distributions’. It is well known that the 
only Gaussian processes that are self-similar and have stationary increments are the fractional 
Brownian motions. Throughout, we let (R/y(t))t e R denote a fractional Brownian motion with 
Hurst index H G (0,1); this is a zero-mean Gaussian process with covariances given by 
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Co <v(B H (t),B H (s)) =-(\t\ + |s|" -|t-s| ), t,s€R. 

Fractional Brownian motions were first introduced in 1940 by Kolmogorov [2(3] and their rel¬ 
evance was first recognized by Mandelbrot and Van Ness [23], who gave them their name. 
Invariance principles for fractional Brownian motions have a long history, since the seminal 
work of Davydov [13] and Taqqu [HI]. As the limiting objects of stochastic models, fractional 
Brownian motions have appeared in various areas, including random walks in random envi¬ 
ronment [3], telecommunication processes [3fl], interacting particle systems 31], and finance 
0 , just to mention a few. 

Recently, Hammond and Sheffield 0 proposed a simple discret model that scales to frac¬ 
tional Brownian motions with H > 1/2. This model, to be described below, can be interpreted 
as a strongly correlated random walk with ±1 jumps. As the simple random walk can be 


Date: April 21, 2015. 


1 





2 


HERMINE BIERME, OLIVIER DURIEU, AND YIZAO WANG 


viewed as the discrete counterpart of the Brownian motion, the correlated random walks pro¬ 
posed in 0 can be viewed as the discrete counterparts of the fractional Brownian motions for 
H > 1/2. In this regime, the fractional Brownian motion is well known to exhibit long-range 
dependence [39]. 

The aim of the present paper is to generalize the Hammond-Sheffield model to dimension 
d > 2 and to study the scaling limits. Based on a natural generalization of the Hammond- 
Sheffield model, we establish invariance principles for a new class of operator-scaling Gauss¬ 
ian random fields. The operator-scaling random fields are generalization of self-similar pro¬ 
cesses (H|) to random fields, proposed by Bierrne et al. Q. Namely, for a matrix E with 
all eigenvalues having real positive parts, the random field (X(t)) teR d is said to be ( E,H)~ 
operator-scaling for some H > 0, if 

(2) ( X(\ E t)) tmd & A H (X(t)) temd for all A > 0, 

where X E := ^ fc>0 (logA ) k E k /k\. In this paper, we focus on the case that E is a d x d 
diagonal matrix with diagonal entries /3i,...,/3d, denoted by E = diag(/?i,..., /3d)- It is 
worth mentioning that a simple generalization of the self-similarity would be to take E being 
the identity matrix in ([21), and the advantage of taking a general diagonal matrix is to be able 
to accommodate anisotropic random fields. Examples of operator-scaling Gaussian random 
fields include fractional Brownian sheets [3] and Levy Brownian sheets [id|. Here, our 
results provide a new class to this family with corresponding invariance principles. We also 
mention that there are other well investigated generalizations of fractional Brownian motions 
to Gaussian random fields, including distribution-valued ones. See for example [3. (25l. [I'D] . 

We now give a brief description of the Hammond-Sheffield model and its generalization 
to high dimensions. Let us start with the one-dimensional model. Let /i be a probability 
distribution with support in {1, 2,...}. Using the sites of Z as vertices, one defines a random 
directed graph Q^ by sampling independently one directed edge on each site. The edge starting 
at the site i G Z will point backward to the site i — Z*, where Z t is a random variable with 
distribution p. Here, g is a probability distribution in form of 


(3) 


= n a L(n), 


where L is a slowly varying function and a G (0,1/2). This choice of a guarantees that 
the graph Q fl has a.s. infinitely many components, each being a tree with infinite vertices. 
Conditioning on one then defines (Xj)j & i such that 

• Xj = Xi if j and i are in the same component of the graph, 

• Xj and Xi are independent otherwise, and 

• marginally each Xi has the distribution (1 — p)5-i + p5\ for some p G (0, 1). 

The partial-sum process S n = Y]” =1 Xi,n > 1, can be interpreted as a correlated random 
walk. Hammond and Sheffield [17l . Theorem 1.1] proved that 


(4) 


f ‘S'lrrtJ — 

\ n“ +1 / 2 L(n) 


C7 (B a+1 / 2 (t)) [ ] 
te[ o,i] 1 J 


as n -MX). Hammond and Sheffield 0 actually established a strong invariance principle for 
this convergence. 


To generalize the Hammond-Sheffield model to high dimensions, we start by constructing 
a random graph Q fl with vertices Z rf . Similarly, at each vertex i G Z d we first sample indepen¬ 
dently a random edge of length Z*, according to a probability distribution (jl, and connect i 
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to i — Z}. The distribution /z has support within {1, 2,... } d , intuitively meaning that all the 
edges are directing towards the southwest when d = 2. Most importantly, the distribution ^ 
is assumed to be in the strict domain of normal attraction of (E , v), denoted by /z £ 'D(E, u), 
for some matrix E = diag(l/ai,..., l/a^) with a* € (0,1), * = 1 and an infinitely 

divisible probability measure v on M+. That is, if (£j)i>i are i.i.d. copies with distribution /z, 
then 

n 

(5) 

i= 1 

This assumption is a natural generalization of (j3j) to high dimensions. We again focus on 
the case that has infinitely many components, which turns out to be exactly the case 
that q(E) := trac e(E) > 2, and given we define {Xj)j £Z d similarly as in dimension one. 
Remark that q(E) > 2 is trivially satisfied for d > 2, due to the restriction on a* € (0,1). See 
Section [2] for detailed descriptions of the measure /z, the random graph Q fl , and the model. 

To study the scaling limit of the partial sums of the random fields over increasing rectangles, 
we introduce 

S n (t) := ^2 x j; n = (ni,...,n d ) £N d ,£ = € [0, l] d 

jeR{n,t) 


with R(n,t ) = rijt=i[0> — 1] H 7L d . Surprisingly, the picture is much more complicated 
in high dimensions. In order to obtain an invariance principle for S n (t ), one cannot simply 
require nrin^-i d n* —> oo as most of the limit theorems for random fields do (see e.g. a i M, 
Ei]). Instead, one needs to investigate 

(6) Sf (t) := X! 

j&R(n E ' l.t) 




with a diagonal matrix E' = diag(/3i,..., /3d)- 

Our main result, Theorem [H reveals the following surprising phenomenon. For different 
E', the limiting random field may not be the same. However, in the special case with E' = cE 
for some c > 0, the dependence structure of the limiting random field is determined by the 
measure v. This case is referred to as the critical regime. For the non-critical regime, one 
can still obtain invariance principles under different normalizations depending on both E and 
E' , although the limiting random field has degenerate dependence structure (either invariant, 
i.e. completely dependent, or independent increments) along at least one direction. Below 
we briefly summarize the phenomenon of critical regime. To the best of our knowledge, the 
existence of such a critical regime has been rarely seen in the literature, except for the recent 
results by Puplinskaite and Surgailis [H, [34| • They investigated invariance principles for a 
different random field model in dimension 2, and referred to the same phenomenon as the 
scaling-transition phenomenon. 


Critical regime: Here we refer to the case of taking E' = E in ([6]). 


Theorem 1. Assume that /z £ D(E,u) for some E = diag(l/ai,..., 1/ad) with on £ 
(0,1), i = 1 ,...,d, and a probability measure v on M^. Assume a\ < 1/2 if d = 1. Let 
if be the characteristic function of v. Then, 


f S E (t)-ES E (t) 


(W 7 (£))te[o,i] d J 

te[o,i] d 
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in the space D([0, l] d ) ; where the limit Gaussian random field (W(i)) tgK d has zero-mean and 
covariance function 



where an explicit expression of a x is given in m below. 

The limit Gaussian random field is easily seen to be (E, H)-operator-scaling with H = 
1 + q(E)/ 2. For this new class of random fields, we study its increments and the Holder 
regularity of the sample paths in Section [5j 

Non-critical regime: For the case E' in Q is not a multiple of E, the situation becomes 
much more subtle. One can still obtain invariance principles with appropriate normalization 
depending on both E and E'. However, in the non-critical regime the limiting random fields no 
longer reflects fully the long-range dependence inherited from Q^. In particular, the covariance 
function of the limiting random field becomes degenerate in certain directions: along these 
directions, the covariance function becomes the one of a fractional Brownian motion with 
either H = 1/2 (the standard Brownian motion, which is memoryless) or H = 1 (the case 
of complete dependence with W(t) = tZ, t > 0 for a common standard Gaussian random 
variable Z). Accordingly, along these directions the increments of the Gaussian random 
fields are independent or translation invariant, respectively. A general invariance principle 
is established in Section U and properties of the limiting random fields are investigated in 
Section [5j Here we only state the invariance principle for d = 2. In the non-critical regime, 
the limit Gaussian random field is a fractional Brownian sheet with Hurst indices Hi and H 2 . 
However, we do not see fractional Brownian sheets in the limit in high dimensions most of 
the time. A complete characterization of when fractional Brownian sheets arise is given in 
Proposition [6] below. 

Theorem 2. Assume d = 2. Let p £ T>(E,u) with E = diag(l/oi, l/af) and set E' = 
diag(l/oi, l/a^) with 01,02 € (0,1), 02 a 2 . Then, depending on the relation between 

01,02 and a' 2 , the following weak convergence may hold: 



in the space _D([0, l] 2 ), where the limit Gaussian random field (W(t )) tgR 2 has zero-mean and 
covariance function in form of 


(7) Cov(IT(£), W(s)) = a 2 X (j 2 Co v(B Hl (H), B Hl (si)) Cov(B H2 (t 2 ), B H2 (s 2 ))- 


Here, (5, a 2 , H\, H 2 and hence {H / (£)} ig [ 0 ,i] ti depend on 01,02 and a 2 . In particular, there 
are four different possibilities as follows: 



The explicit expressions of a 2 in these cases can be found in the proof of Theoroem 0 in 
Section [5j 
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The main result of the paper, Theorem [5l is a unified version of invariance principles for 
general rfeN,l? = diag(l/ai,..., 1 /ad) and arbitrary E 1 , from which both Theorems Q] and [2] 
follow as immediate corollaries. Theorem [5] also provides a general principle to determine 
the correct normalization order, the limit covariance function, and hence the directions of 
degenerate dependence. We have just seen that in dimension 2 there are already 4 different 
non-critical regimes. For general d > 3, the situation becomes more complicated. 

The core of the proofs is an application of the martingale central limit theorem, thanks to 
the key observation that the random field of interest can be represented as a linear random 
field in form of 

(8) x i =Y J qjX*i- j ,ie z d , 


of which the innovations {Xj)j eZ d are multiparameter martingale differences. Hammond and 
Sheffield [17|] also made essential use of the martingale central limit theorem, although the 
representation as a linear process as in © was not explicit. This representation plays a key 
role in our proofs, as from there when verifying conditions in the martingale central limit the¬ 
orem, thanks to the structure of the linear process, we can deal with the coefficients qj and 
innovations separately. This framework, or more generally the martingale approximation 
method, has been carried out successfully in dimension one to establish invariance princi¬ 
ples for fractional Brownian motions for general stationary processes jT^]. To extend this 
framework to high dimensions, a notorious difficulty is to find a convenient multiparameter 
martingale to work with. It is well known that the martingale approximation method applied 
to stationary random fields is not as powerful as to stationary sequences, as pointed out long 
time ago by Bolthausen [§|. Fortunately, our specific model can be represented exactly as a 
simple linear random held with martingale-difference innovations as in ([8]) . 

Once the representation of linear random fields in (JS]) is established, the main work lies 
in the computation of the limit of the covariance functions. This step is heavily based on 
the analysis of Fourier transforms of the linear coefficients (gOieZ'D the asymptotic property 
of which is essentially determined by v. Analyzing the Fourier transforms is a standard 
tool to compute the covariance functions for stationary linear random fields, see for example 
[23:, |33j, lH]. To complete the invariance principle, the tightness is established. At last, to 
develop the sample-path properties we apply recent results in Bierme and Lacaux jjj. 

The rest of the paper is organized as follows. In Section [2] we describe in details the random- 
held model. Section [3] provides a general central limit theorem that serves our purpose. 
Section [4] establishes a general invariance principle that applies to both critical and non- 
critical regimes. Some properties of the limit random helds are provided in Section [5j 


Acknowledgement. The third author would like to thank the hospitality of Laboratoire 
de Mathematiques et Physique Theorique, UMR-CNRS 7350, Tours, France, during his visit 
from April to July in 2014, when the main part of this project was accomplished. The third 
author’s research was partially supported by NSA grant H98230-14-1-0318. 


2. The model 

In this section, we will give a detailed description of our random held model {Xi} ieZ d, 
of which the dependence structure is determined by an underlying random graph Q^. The 
asymptotic properties of the random graph is determined by a probability measure /1 on 
{1,2,... } d , which is assumed to be in the strict domain of normal attraction of an Fl-operator 
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stable measure v on M^. Some simple properties of the model will be derived. In particular, 
we show that the random field of interest can be represented as a linear random field, of which 
the innovations are stationary multiparameter martingale differences. 

Throughout the paper we use the following usual notations. Let d > 1 be an integer. 
On M. d , we consider the partial order (also denoted by <) defined by t < s if tj < Sj for 
all j = 1,..., d, where t = (t\,... ,td) and s = (s i,..., sfi). In the same way, we use the 
notations >, <, >. We write t s as soon as tj > Sj for at least one j = 1 ,,d, and in 
the same way, we use We denote by [£, s] the set [fi, si] x ■ ■ ■ x [f^, s<i] and we write 

|i|oo for max{|tj|, j = 1 ,...,d}, and |t|i for Xq=i l^'l- Furthermore, write N = {0,1,...} 
and N* = {1,2,... }. 


2.1. The random graph. 

Let gbea probability measure on N* such that the additive group generated by the support 
of p, is all Z d (we say that fj, is aperiodic). On Z d , we consider the random directed graph Gp, 
associated to /r, defined as follows: 

• Let (Z n ) neZ d be i.i.d. random variables with distribution \i. 

• For each n£Z{ let e„ be the outward edge from nton- Z n . 

• G/i, is the graph with all sites of Z rf as vertices and random directed edges {e„, n £ Z d }. 

The graph G^ is then composed of (possibly) several disconnected components and each 
component is a tree. The upcoming Proposition [T] shows that, almost surely, the number of 
components of Gp, is one or is infinite. 

We first introduce the following notations. For n € Z d , we denote by A n the ancestral line 
of n, that is the set of all elements k £ Z d for which there exists a directed connection from 
n to k (taking the orientations of the edges into account). Note that, in distribution, A n 
can be described by the range of the random walk (n — Sk)k >o where (Sk)k >o is the random 
walk starting at 0 with step distribution / j . In particular, since fi is supported by N d , any 
element k in A n satisfies k < n. Observe that the condition that the support of generates 
the group Z rf is equivalent to the fact that P(A n 0 A m ^ 0) > 0 for all n, m £ Z rf . 

For n £ Z d , we set q n = P(0 £ A n ). We clearly have q n = 0 as soon as 0 ^ n, except for 
go = 1. Further, since each edge is generated independently at each site, for any n, k £ Z d , 

P(fc £ A n ) = q n -k- 

Proposition 1. IfJ2keN d converges, then Gp has almost surely infinitely many components 
whereas if Qk diverges, then Gp has almost surely only one component. 

We start by proving the following lemma. 

Lemma 1. (i) converges then for all n £ Z d , 


''f ] QkQk+r 


PO4 o nA^0)= J] 9* 

\fcg N d 

(ii) If Ik diverges then P(Ao 0 A n 0) = 1 for all n £ Z d . 

Proof. The proof follows an idea developed in fl7l . Lemma 3.1] for the dimension 1. Let G', 
be an independent copy of Gp- We denote by A' n the ancestral line of n with respect to Gp- 
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On one hand, one has 

E|ALo O A' n \ = P(& £ Ao)P(fc £ A n ) = QkQk+n- 

k&z d feez d 

On the other hand, 

E|A 0 O A'J = P(A 0 0 A n ? 0)E|A o 0 A' 0 | = P(A 0 0 A n + 0) ^ q 2 k 

fee N d 

and thus jG|) follows. 

If ^fceN d ( lk = °°) then E|Ao O A' 0 \ = oo. But E|Ao O A' 0 \ can also be computed as 

(9) E|A 0 O A' | = ]T F(\Ao n 4,| > fc) = £ E(^o 0 A' 0 jt { 0 }) fc = 1 my 

k> 0 fc>0 ' 0 oT\ t) 

Thus EIAoOAqI = oo if and only if P(AoOAq / {0}) = 1, and in this situation |AonA' 0 | = oo 
almost surely. Now, since the group generated by the support of // covers Z d , we know that, 
for all n £ Z rf , there exists fco £ Z rf such that 

P(fc 0 ^ A 0 and k 0 — n £ A' 0 ) = P(fc 0 G A 0 0 A' n ) > 0. 

But, since |Ao 0 A' 0 \ = oo a.s., we infer that | A ko 0 A ko -n\ = oo also a.s., and thus 

P(Ao ni n /0) = P(A fco n A ko _ n + 0) > P(|A fco n A fco „ n | = oo) = 1, 

which proves ImI) . 

Proof of Proposition [7J If C := ]C fcgN d < oo, from Lemma Q] |zj), we get 


□ 


P(A 0 ni n /0) = cr 1 ^ g fc%+n < C -1 


ke%j d 


<k€.N d , k> —n 


E 

k>n 

oo, and we can build a 


which goes to 0 as (nloo — y oo. Thus, P(Ao H A n / 0) — > 0 as |n|o< 
sequence {n k ) k& ^ C Z d , iteratively, such that for each k £ N, 

P (Ai fc n (U j =0 A nj ) / 0^ < p- 

By the Borel-Cantelli lemma, we see that, almost surely, the ancestral lines A nk , k > 1, are 
disjoint from each other. This proves the first part of the proposition. 

The second part of the proposition is clear from Lemma [T] / Tm]) . □ 


2.2. The measure. 

From now on, we always consider an aperiodic probability measure p on N* such that 
P £ D(E,v), as defined in ([5]), where v is a full probability measure on and E = 
diag(l/ai,..., l/arf). Since the distribution of each coordinate is in the strict domain of 
normal attraction of a positive stable laws and since positive a-stable laws only exist for 
a £ (0,1), it necessarily follows that an £ (0,1) for all* = 1,..., d. 

In this case, p is also said to have non-standard multivariate regular variation with exponent 
E. Equivalently, p is (non-standard) multivariate regularly varying with exponent measure (f, 
with cj) being the same Levy measure in the triplet representation of v [29|, Corollary 8.2.11]. 
That is, for some constant c > 0, 

(10) lim np{n E A) = af>(A) for all A £ £>(M d ) bounded away from 0 and (f>(dA) = 0. 

n—>oo 
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Equivalently, this means that np(n E -) converges vaguely to ccf>, in the space of Radon measures 
on \ { 0 } equipped with the vague topology, under which sets in F(R d ) bounded away from 
{0} are relatively compact. Most of the applications in the literature of multivariate regular 
variation, however, focus on the case that aq = ■■■ = ad- In this case, cnD is referred to as 
multivariate regular variation in the literature. Standard references on (standard) multivariate 
regular variation includes [32, 38]. References on non-standard multivariate regular variation 
include [35j], (38, Chapter 6]. See also some recent development in [36]. Some examples are 
given at the end of the subsection. 

We denote by P the Fourier transform of the measure /i, that is 

P(t) = Y ^{{k}y t k , t G R d . 
fee N d 

Note that the assumption that the additive group generated by the support of p, is all 7L d is 
equivalent to: 


P(t) = 1 if and only if the coordinates of t belong to 27 tZ, 
see for example Spitzer (13, p.76]. 

Let Qn be the random graph associated to p as defined in Section 12.11 The asymptotic 
behavior of {qk}keN d pl a Y a Ley r °fo our analysis. It is essentially determined by the 
measure p G T>(E,v). We denote by Q the Fourier series with coefficients qk = P(0 € Ak), 
that is 

Q(t) = Y Qke it k ■ 
fee 


Using that qk 
the relation 


T({j})qk-j for k > 0, we see that both Fourier series are linked by 


Q(t) 


l 

i -p(ty 


From Lemma [TJ we see that 


P(^o n A n + 0) 


c n (\Q\ 2 ) 

co(|Q| 2 )’ 


where Ck(\Q\ 2 ) denotes the Fourier coefficient of index k of \Q\ 2 = QQ. This relation explains 
why the Fourier series Q plays a crucial role in the study of the random graph. 

We denote by ijj(t) = J Rd e lt x di'(x) the characteristic function of the full F-operator stable 

measure v. Note that it follows from Q that the log-characteristic function logi/> is then an 
F-homogeneous function, that is 


for all f > 0 and x G M rf , log i))(t E x) = t log ijj(x). 


Further, logV’(O) = 0 and for all x / 0, | log^(x)| > 0. 

The two following lemmas are key results concerning the behavior of Q at 0. 


Lemma 2. Let fi G D(F, v) be as described above and i/j the characteristic function of v. 
Then 

iq(*)I = ii - p(x) r 1 = x e 

where g is continuous and positive with g( 0) = 1. 
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Proof. Let us use a change of variables in polar coordinates. As in 
a new norm on M rf , related to the matrix E, by 

1 


( 11 ) 


1*11 E = 


f 


E , 
r x 


-dr, 


Chapter 6], we define 


where here | • | denotes the Euclidean norm. The unit ball Se = {x £ | j|*||e = 1} 

associated to this norm is a compact set of \ { 0 } and every vector in \ { 0 } can be 
uniquely written as r E 0 with r > 0 and 0 £ Se , since for any x ^ 0, the map t i-A ||f E a:||£ is 
strictly increasing on (0, oo). 

Since p £ T>(E, v), we have 

P(n~ E 0) n —>• ip(0), as n — » oo, uniformly in 6 £ Se, 


from which we infer that 


t\ogP(t~ E 0) -A- \ogif{6), as t —>• oo, uniformly in 6 £ Se, 
see [H, p. 159]. Using that log(l + x) ~ x as x —> 0 and that P is continuous at 0, we obtain 
t.(P(t~ E d) — !)—>• logip(G), as t —> oo, uniformly in 6 £ Se- 


Thus, for all e > 0, there exists T > 0 such that for all t > T, 


log ip(t E 0)\ 


| P(t~ E 0) - 1| 


- 1 


I log ip(9)\ 


t\P{t~ E 0) - 1| 


- 1 


< e, uniformly in 6 £ Se- 


Now, set g(-) = \ log — 1) 1 |. The function g is clearly continuous and positive on 


[—7r, ir} d \ {0}. Set 5 = inf o & s E 11^ h ^\\E > 0- Then for all x such that ||*||_e < 5, x = t 0 E 6o 
with <?o £ Se and to > T and thus 

\g(x) - 1| = \g{tQ E e 0 ) - 1| < e. 

Thus g is continuous at 0 and (?(0) = 1. □ 


We are thus interested by the function x i —> log ip{x), which is a continuous E-homogeneous 
function that only vanishes at 0. Recall that q(E) = trace(E'). 

Lemma 3. If <f> : —> M is a continuous E-homogeneous function that only vanishes at 0, 

then for any p > 0, x t —> \cf(x)\~ p is locally integrable in if and only if q(E) > p. 


Proof. There exists a unique finite Radon measure oe on Se which allows the change of 
variable 


t- r+oo 

/ = / 

jR d Jo 


t Se 


f{r E d)A E ^- l da E {e)dr , 


for all / £ L 1 (W i ) (see [6j, Proposition 2.3). Thus, using the E-homogeneity of cf, one has 
f \(j)(x)\~ p dx = / / r q ^ E ^ l \(j){r E 6)\~ p d(TE{0)dr 

J {||a31|£7<1} Jo JSe 


f 1 r q(E)-i- Pd r f \^q)\- p da E {6). 
Jo Js E 


The second integral is finite because \(j)\ is continuous and positive on the compact set Se, 
and the first integral is finite if and only if q(E) > p. □ 


As a first consequence, we get the following proposition. 
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Proposition 2. Let p G T>(E,v). The random graph Qf has almost surely infinitely many 
components if and only if q(E) > 2. 


Note that, when d = 1, the condition q{E) > 2 becomes a\ < which corresponds to 
the condition assumed in [17]]. When d > 2, since a* G (0,1) for all i = 1,... ,d, then the 


conditon q{E) > 2 is always satisfied. 


Proof. As a consequence of Lemma [21 using Parseval identity, we get 


J2 

fee 


1 

(27 l) d 


f \Q{x)\ 2 dx 

J [—7T,7r] d 



| 5 , (*)| 2 |logi/>(ai)| 2 dx. 


Since g is bounded and bounded away from 0 on any compact set, we see that q\ < +oo 

if and only if a; i — >■ | log ^(*)|~ 2 is integrable on [— 7r,7r] d . The function x i —> \ogfi(x) being 
E-homogeneous, by Lemma [21 it is the case if and only if q(E) > 2 and the result follows 
from Proposition [H □ 


To conclude the section, we give few examples of possible probability measure p G T>(E, u). 


Example 1 (Product measure). Let p be the product measure pi <S> ■ ■ ■ <S> Pd, where each pi 
is a regularly varying measure on N* with index G (0,1) such that 

Pi([n, oo)) ~ Cin~ ai , 

for some <7 > 0. Then, each pi belongs to the strict domain of normal attraction (with 
normalization n -1 / ai ) of a positive Oj-stable law z^, see [3, Theorem 8.3.1]. Positive a-stable 
laws only exist for a G (0,1), and then, their characteristic functions are given by 

pit) = exp | — 7 |t|“ (l - isgn(t) tan (^a))} , 

for some 7 > 0. See [3, Theorem 8.3.2]. In this situation, the measure p belongs to the strict 
domain of normal attraction of the measure v = 17 (g> • • • (g> Vd which is a full E-operator stable 
distribution, with E = diag(l/«i,..., 1 /ad). The characteristic function if of v is such that 

d 

log 'fiix) = 5 ^ 7 j\xjp (l - zsgn(xj) tan , 

j= 1 


for some 7 j > 0 . 


Example 2 (Standard multivariate regular variation). For the standard multivariate regular 
variation, that is when ot\ = ■ ■ ■ = ad = at, many examples have been known from the studies 
of heavy-tailed random vector X = {X\, ..., Xfi) G M rf , in the literature of heavy-tailed time 
series. An extensively investigated condition for multivariate regular variation is 


( 12 ) 


P(|X| > ux,X/\X\ G •) 
P(|X| > u) 


=7 x a cr(-) as u —> 00 , for all x > 0, 


for | • | a norm on and a a probability measure on B{S) for S = {x G : |x| = 1}. See 
for example |l]. It is known that (HOD implies (I12D (see e.g. [2d3 . Theorem 1.15]). 

The measure a is often referred to as the spectral measure, which captures the dependence 
of extremes. For example, the case that a concentrates on the d-axis with equal mass means 
that, in view of (I12D . the extremes of the stationary processes are asymptotically independent. 
For more theory and examples on spectral measures reflecting asymptotic dependence of the 


extremes, we refer to [37|, Chapter 5]. 
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Example 3 (Polar coordinate). A standard procedure to obtain non-standard regularly 
varying random vectors is via the representation using polar coordinate. We use the 
norm || • \\e introduced in (flTT) to identify \ {0} with (0,oo) x Se for the unit ball 
Se = {x £ | \\x\\e = 1} such that every vector in \ {0} can be uniquely written 

as r E 0 with r > 0 and 6 £ Se- By [29|, Theorem 6.1.7], in case of ([5]) (equivalently (flOlP . 0 
can be taken to have the polar coordinate representation 

r°° r dt 

(13) 0(A) = / / t {tE0eA} a(dO)-2, 

J 0 JSe “ 

for some finite Borel measure a on Se- In our case, since p has support contained in N*, 0 
is a measure on K(|_, and a is a finite measure on S^ = Se D M+. Identifying \ {0} with 
(0, oo) x Sfi, to obtain a multivariate regular varying measure as in (ffoD, it suffices to show 

(14) /r((r, °°) xT)~ cr _1 (j(r) as r —> oo, for all T € B(S^). 

This follows from a standard argument showing that {(r, oo) x V } r>0 rg g( 5 +) are a convergence 
determining class. 

A standard procedure to construct a random vector of which the distribution p satisfies (1141) 
is the following. Let R be a non-negative random variable with P (R > r) ~ ca(S'^)r~ 1 as 
r —> oo. Let © be a random element in 5^ with probability a/a{S g). Assume that R and 
0 are independent. Then, R E & is regularly varying in M'j. in the sense of (jl4[) . Indeed, 

P (R e & £ (r, oo) x T) = P(i? > r, 0 £ T) ~ cr _1 ir(r) as r —>■ oo. 

The so-obtained distributions can then be modified to become distributions on N* with the 
same regular-variation property. We omit the details. 


Remark 1. For our main results to hold, we do not impose any assumption on the spectral 
measures in Examples [2] and [3] The only assumption is the non-standard multivariate regular 
variation with indices a\,... ,aa £ (0,1), and aq < 1/2 when d = 1. 


2.3. The random field. 

We now associate a random field to the random graph Q^. Assume that /U £ 

~D{E, v) as in the preceding section, with the diagonal matrix E satisfying q(E) > 2, and let 
p £ (0,1). We proceed as follows: 

First, generate the random directed graph Q^ as described in previous sections, which has 
almost surely infinitely many connected components in this situation. Let {Ci \ i > 1} denote 
the collection of disjoint components and associate to each component Cj a random variable 
£i such that i are i.i.d. with distribution given by P(ej = 1) = p and P(ej = —1) = 1 —p. 
Finally, for all j £ Z d , set Xj = e, where i is such that j £ C*. This construction implies that 
Xj = Xj~ as soon as j and k belong to the same component of Q fi , and they are independent 
otherwise. 


Remark 2. As pointed out already in 0 , the one-dimensional model is an example of the 
so-called chains with complete connections, which has a long history with different names; 
see flR 16 ] for more references. In the same spirit, our model is an example of partially 
ordered models recently introduced by [l3[ , an extension of chains with complete connections 
to random fields. 
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For all n £ N d , we introduce the partial sum 



Our aim is to establish a functional central limit theorem (invariance principle) for the 
partial sums S n (with centering and appropriate normalization) when n goes to infinity with 
a specific relative speed in each direction. We will distinguish different regimes. We first 
show, in this section, that (Xj)j e %d. can be seen as a linear random held with martingale 
differences innovations, and thus, S n is a partial sum of a linear random held. 

For all j £ 7L d , we define the cr-helds = cr {X^ | k < j} and Wj = a{Xk | k ^ j}. Note 
that, for j < n, the value of X n conditioned on Gj is obtained by sampling the ancestral 
line A n and taking the value of X & where k is the hrst site of the ancestral line A n which is 
strictly smaller than j. We denote 


x; = Xj - E (Xj I <Tj) = Xj - E(Xj I Wj). 


(15) 


The equality E(Aj | Wj) = E (Xj \ aj) comes from the fact that starting from j, the next 
site in the ancestral line Aj is necessarily strictly smaller than j . Then for all j £ E d , 
E (Xj | Wj) = 0 and Xj is measurable with respect to ~Vj+ eq for all q = 1,..., d, where e q is 
the (/-th canonical unit vector of M d . In particular, the random variables Xj are orthogonal 
to each other, that is, E(XjX^) = 0 as soon as j / k. 

Lemma 4. In the above setting, 



Proof. Let Zo be the random variable with distribution /r that gives the first ancestor of 0. 
We have X 0 = J2k>o t {z 0 =k} x -k and E(A 0 |cr 0 ) = J2k>oPk x -k, where p k = p({k}) for all 
k > 0. Therefore, 



(16) 


^2 “ Pk)(^{Z 0 =(.} - Pi))^{X^ k X-l). 


k> 0 £>0 


But, 

(17) E(X_ k X_ £ ) = P (A_ k n A_ £ ± 0)E(^2) + P(Al_ fc n A_ £ = 0)E(A' o ) 2 . 


and 


(18) 


E ((l{Z 0 =fc} ~Pk)(l{Z 0 =£} - Pi)) = l{fc=£}Pfc ~PkP£- 
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Combining (fl6|) . (fT71) . and (fT8l) . we get 
E(X^ 2 ) 

= E(X 2 ) L-EE PkP(F(A_ k n / 0) -EE PkPt^(A_k n A_£ — 0)E(Xo)" 

V k>0£>0 / k>0l>0 

= 0E(X 2 )-E(X O ) 2 ) EE PkPfF(A_ k n = 0) 

fc> 0£>0 

= Var(X 0 )P(A 0 n X = {0}), 

where A' 0 is an independent copy of Aq. Finally, as we saw in ([9]) in the proof of Lemma 
= H A' 0 \ = E(^4o H A' 0 = {0})” 1 and the proof is complete. 

Now, for all j E Z d , we introduce 

A j(X)= J2 (-l) d " |£|l IE(A | a j+e ), 

e€{0,l} d 

where |e|i = £\ + ... + 

Remark that, since E (Xj \ Vj+e) = ^(Xj \ aj ) for all e E {0, l} rf with the exception of 
e = 1 for which E (Xj \ crj+i) = Xj, we have 

(19) A j{Xj) = Xj - E {Xj \ aj) = X*. 

More generally, we have the following lemma. 

Lemma 5. For all j, k E Z d ; 

A j(X k ) = q k _jXj , 

which vanishes when k f> j. 

Proof. The result is clear when k = j (see (119(1 . In the case k < j, k ^ j, we easily see that 
A j(X k ) = 0. 

Now, assume k ^ j. By linearity, we have 

A j{X k ) = Aj(X k l{j £ A k y) + Aj(X k l.{j£A k }). 

Using first that X k t{j^A k } = Xj^{j£ Ak y, and then that {j E A k } is independent of (Tj+i, we 
obtain 

Xj(X k l{j eAh }) = Aj (Xj )P(ji E A k ) = q k _jXj. 

Denote by a(k,j) the hrst element of the ancestral line A k that is < j and remark that 
a(k,j) is independent of crj+i- Then, 

Xj(X k t{j£ Ak }) = Xj(X k l{ a(k j) = i}) = A j{Xe)F(a(k,j) = £). 

But, A j(Xi) = 0 for all i < j, t ^ j, and we finally have 

Aj( x k t{j<£ Ak} ) = 0 , 

which completes the proof. □ 

Lemma 6. For all k E Z rf , the series A j(X k ) converges in L 2 and 

Xk- E(X fe ) = 22 Ai(X fc ). 


3 


□ B 
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Proof. First, remark that by stationarity we may only consider the case where k = 0. The 
sum in the statement can be write as Jfj eN d A_j(Xo) since the other terms vanish. We 
denote by nl the vector (n,..., n) where n £ N. By Lemma 01 we have 

E ( { E ) = E(-V 2 ) ( E 

\ \ie[o,ni] / ) \je[o,ni] / 

and the right hand side converges to Var(Ao) as n —>• oo thanks to Lemma 01 Now, by 
construction, the random variables X)je[0nil ^—j( x o) and Xo — X)je[o nil ^-j(X o) are or- 
thogonal. To see this last fact, note that for all l < 0 and j < 0, E (E(X 0 | ai) | <Jj) = 
E (Ao | cr min {7 j}) , where the minimum is taken on each coordinate. Thus, we get 


E 


Xo-E(Xo)- Y = Var(A 0 )-E ]T 


ie[o.ni] 


yje[0,nl] 



as n —> oo. 


□ 


From Lemma [6] and Lemma 01 we get that ( Xj — E (Xj))j eZ d is the linear random field 
given by the innovations (Xj)j €Z d and the filter (qj)j & z d - That is, for all k £ Z d , 

X k - E(X fc ) = Y <lk-jX*. 

Hence, we proved the following proposition. 

Proposition 3. For all n £ 

S n - E(5„) = Y <3 X *j- 

j£Z d 

where b n j = I] fce [ 0 , n _i] Qk-j■ Further, for any n £ N d , b n = (b n j) j& d belongs to l 2 (Z d ), 
that is \\b n \\ 2 := YljeZ d ^i,j < oo. 


3. A CENTRAL LIMIT THEOREM 

We still assume fi £ D(E,u), where u is a full El-operator stable law on with E = 
diag(l/ai,..., 1/ad), with cij £ (0,1) and a\ £ (0,1/2) if d = 1. The random field (Xj)j £Z d 
is the random field defined in Section 12.31 In view of Proposition 01 we want to establish 
central limit theorems for the sequences of L 2 random variables 

Y, °n, X h 71 ^ 1 

jez d 

with general coefficients c n = {c n ,j)j & z d £ ( 2 (1‘ d ). Recall the definition of AT in (1X5]) . It turns 
out that a simple assumption on c n for a central limit theorem is given by 

(20) lim sup = 0. 

n ^°°jez d IMI 

The aim of this section is to prove the following central limit theorem. 
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Theorem 3. Let c n 


where 

( 21 ) 


( c n,j)j£Z d be a sequence in( 2 (h d ) satisfying ([20]) . Then 



Y,c nJ X^Af(0,a 2 x ) 

j£Z d 


as n —>• oo, 


4 := Var(X*) 


Var(Xp) 

Tk 


Proof. Recall that we write aj = a{X /, k < j} and aj = a{Xj~ \ k ^ j}, and we already 
have seen for all j G Z rf , 

I a j ) = I a j)- 

We now consider the c-fields Tj = a{Xj, \ k -< j}, where -<; denotes the lexicographical order 
on 7i d . We have aj C Tj C aj for all j G T, d and thus, for all j G Z rf , we also have 


E (Xj | Tj) = E(Xj | aj). 


In general, if {Ti}i & id is a filtration such that Ti C Tj if * -< j. for all i,j G Z d , we say 
that integrable random variables (£i)i g z d are martingale differences with respect to {-T*}i e z d 
if 


fi G Ti +ed and E(£j | Ti) = 0 for all i G Z rf , 
where is the d-th vector of the canonical basis of 

Thus, by definition (see (fl5]l ). the random field {Xj)j £Z d is composed of martingale differ¬ 
ences with respect to the filtration {Tj}j & %d defined above. As a consequence we will be able 
to use the following theorem which is an obvious adaptation of a theorem of McLeish [2£*] for 
triangular array of Z-indexed martingale differences. 


Theorem 4 (McLeish [28)). Let (f,n,j) n eN,jez d be a collection of random variables satisfying 
Y!j£Z d fin ) j £ L 2 for all n G N. Assume that for each n G N, ( [f,n,j)jez d cire martingale 
differences with respect to a filtration {T n j } 3 - 6Z d in the lexicographical order. If 

(i) limj^oo max^g^d \f n ,j\ = 0 in probability, 

(ii) sup neN E (max ieZd < oo, 

(Hi) lim n _>.oo 52jezd(,n,j = c 2 > 0 in probability, 
then 

=> A/"(0, a 2 ) as n —>• oo. 

j£zL d 


Proof. Let us explain how one can derive this theorem from Theorem 2.3 in [28] 
stated for finite sets of random variables at each n. First, since Y)jez d £n,j A 2 , 


find a sequence of finite rectangles in Z rf such that J2jez d \r n £n,j converges to 0 
n —> oo. Thus, the conclusion of Theorem 0] holds as soon as 


which is 
one can 
in L 2 as 


=> a 2 ) as n —> oo. 

jer n 


Furthermore, for each n, using the lexicographical order on the finite set r n , one can re-index 
the random variables (£ n ,j)je r„ and the er-fields {T n j}j £ r n in order to fit with the statement 


of 

of 


281 . Theorem 2.3]. Now, it suffices to observe that conditions 
28, Theorem 2.3]. 




and 


imply those 

□ 
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Theorem [3] will be established by application of Theorem 2] to 

Znj ■= and F nJ := Fj = a{X k \ k -< j}. 

11 11 

Note that \Xj\ < 2 and thus the conditions 0 and |m]) can be reduced to a condition on the 
coefficients c n j. Indeed, @) and (@) are satisfied as soon as (E0|) holds. Condition ifm)l can 
be derived from the following lemma. 

Lemma 7. Let c n = {c n j)j &Z d. be a sequence in l 2 (Z d ) such that (1201) holds. Then, 

lim EE V c 2 n ,Xf = E^ 2 ) in L 2 , 

rwoo \\c„ 2 t—* ’ 3 3 

" " jez d 


Proof. We start by showing that 

(22) Cov(Xf 2 , Xj 2 ) -A- 0, as | i — j\ ao —> oo. 

Observe that X*. = Xj - Y,t>oPt X j~e and let X j,k = x j ~ X^e{i ,...,k} d Pi x j-i- For an y 
j £ Z d , using that \ x *\ < 2, we get 


\Xf - X*A\ < 


V* V* 

\ X j - x j,k I 


= 4 


e pfXj-i 


^G[l,oo) d \[l,fc] <i 


Thus, since \Xj \ = 1 for all j £ Z d , 

(23) sup \X* 2 - X*\\ < 4/i ([1, oo) d \ [1, A:] d ) a.s., for all k > 0. 

Now, introduce 


Ri i k 


u * 


n 


u 


Am = 


\££i— [0,fc] d 


ym£j — [0,k] d 


We have 

nRi,j, k )< E E A m ^%). 

£.£i—[0,k] d m£j—[0,k] d 

But, from Lemma [T] (E), we see that E(A^ O A m ^ 0) —> 0 as \i — m 1^ —>• oo and thus, for 
any k > 1, 


(24) ¥(Rl j k ) —> 0, as |i — j\oo °°- 

Fix e > 0 and, using (1231) . let k £ N be such that supj eZ d |X| 2 — Xj 2 k \ < e. From (1241) . for 
\i — j loo large enough, we have E(Rjj k ) < £ and we obtain 

E{Xfxf) = n x ll x i%) + 0(e) = E(X* 2 k X* 2 k I R iJik ) + 0(e) 

= E(X? 2 | Rij'k) E(X* 2 fe | R. tJ . k ) + 0(e) = E(Xf 2 )E(X; 2 k ) + O(e) 

= E(X? 2 )E(Xf) + 0(e). 


This proves (021) . 
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To prove the lemma, fix e > 0 and let K > 0 be such that | Cov(A? 2 , Xf 2 )\ < e as soon as 
|i — j’Icxd > K. One has, 

n 2 

1 


E 


J2 c 2 nJ x f -E(X* 0 2 ) 


„2 


c 2 ■ 


EZ^E 


j£Z d 


< 


\\r II2 Z_^ 


C 2 ■ 
n.z 


Crr Con 

■' n11 |i-i|oo<AT 11 n| 

r2 

'n.k 


|Cov(Xf,Xf)|+ £ j;^ 


c 2 ■ 

71,2 


|*-i|oo>A: 


< sup ]T \Cov(X* 0 2 ,X?)\+e, 


kGZ d II c n 


|*—0|oo<A' 


and the first term of the right hand side goes to 0 as n —> oo because | Cov(Aq 2 , A? 2 )| is 
bounded and sup fcgZ d c 2 k = o(||c n || 2 ) by assumption. □ 

Thus the conditions (Gj), |m]), and IHM) are satisfied with a\ = E(Xq 2 ) = Var(Xg) = 
(SfceN d Var(Ao) by LemmalU Then Theorem |4] applies and using Lemma0]we complete 
the proof of Theorem [3] □ 

The following lemma gives another useful condition on the coefficients (cn,j)jez d f° r The¬ 
orem [3] 

Lemma 8. If (c n j)j £Z d is a sequence in I 2 {fL d ) that satisfies, for all q = 1,..., d, 

1 


(25) 


lim 

n->oo II Cn 


E \ C l. 


3 c n,j+e q 


= 0 , 


j ez d 


where e q is the q-th vector of the canonical basis of M d ; then m holds. 


Proof. We use an idea of [32]. Assume that (12U1) does not hold. Then, there exist e > 0, a 
sequence {rik)k >l such that -» oo as k —> oo, and a sequence ( jk)k>i such that c nk ,j h > 
e||c n J| for all k G N. Choose M > 0 such that M d e 2 > 1. One has, for all fcsN, 

|| g || 2 > \ ' M d c 2 • _ \ . . _ g 2 | 

II C n k ,j k +j £- 1V1 c n k ,j k \ c n k ,j k +] c n k ,j k \- 


Hence, 

(26) 


je[o,M-i] d je[o,M-i] d 

(mv - ljikji 2 < 23 

j'E[ 0,M—l] d 


2 _ 2 
C nkJk+3 Cn k,3k I' 


But, if j G [0, M — l] d , then 




I 2 — 2 I < 

I C n k Jk ° n kJk+3' ~ 


E i 2 _ 2 


2=1 
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where A = (A 0 , Ai,..., A^) is any path from A 0 = j k to A^ = j k + j, with | A,; — A,; +1 |! = 1 . 
Since j € [0, M — l] rf , we can always choose the path A of length i = ^(A) smaller than dM. 
Thus, we get 


I C nkJk C nk,jk+j\ — dM 


sup sup \c nk)k -c nk)k+eq \ 

q=l,...,d 


d 


— XT \° n k, k c n k ,k+e q 

q =l kez d 


Together with (PItJD . this contradicts (1251) . 


□ 


Remark 3. Using Cauchy-Schwarz inequality, we also see that the condition 

(27) lim ^ V (c n j - c nJ+eq ) 2 = 0, for all q = 1,..., d, 

™ IMI j&d 

implies (1251) and thus by Lemma| 8 l implies (1201) . This last observation leads to an improvement 
in Theorem 3.1 in Bierme and Durieu [3]. The conditions (i) and (ii) of this theorem are 
equivalent to our conditions (l20l) and (l27l) . respectively. Thus, the condition (i) in (3, Theorem 
3.1] is unnecessarily. 


4. An INVARIANCE PRINCIPLE 


The aim of the section is to establish a general invariance principle for partial sums of 
the random field (Xj)j e7i d defined in Section [2j Recall that ( Xj)j £Z d are associated to the 
random graph Q^, with fi € T>(U, v). We consider partial sums on finite rectangular subsets 
of h d . As we will see, the growth of the rectangles will be determinant in the invariance 
principle and different limit random fields appear at different regimes. For the general case, 
consider a matrix E' = diag(l/a / 1 , .... l/c^) with a[ > 0, i = 1,... , d and the partial-sum 
process 

Sj?(t) = ^2 X E n > 1 and t = (ti,... ,t d ) <E [0, l] d . 

j£[0,n B 't-l] 

The result will depend on both E' and E. 

We introduce several parameters. First, for all k = 1,..., d, set := ak/a' k , and consider 


(28) 


7 o = 7 o (E,E') 


min 


7 e Pd} 


k-7>Pk k-l>Pk ) 


Note that 70 is well defined by the assumption q(E) > 2, and is completely determined by E 
and E’. Given 70 > 0, define the sets 


/< :={/;€ {1,... ,d} | 70 < Pk}, 
1= := {k € { 1 ,... ,d} | 70 = p k }, 
/> := {k £ { 1 ,... ,d} | 70 > p k }. 


This gives a partition of {1, ... ,d}. We also write /< := I< U I = and /> := I = U />. The 
sets I> and /< consist of the directions in which the limit random held exhibit degenerate 
dependence structure. Remark that by construction, 


|/=| > 1 and |i>| < 1. 
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According to these subsets of {1,... , d}, we consider subspaces of given by 

'H < := {x G | Xk = 0 for k ^ /<}, 

and similarly H = , Let vr<, 7 r=, 7 r>, 7 r<, and 7 r> denote orthogonal projections to 

the corresponding subspaces, and let A<,A=,A>,A<, and A> denote the Lebesgue measures 
on the corresponding subspaces. For 7 r of any proceeding projection, ttE is a linear operator 
on M d ; accordingly there is a corresponding diagonal matrix, which we also denote by -kE 
with a little abuse of notations. 

Next, we define another diagonal matrix E" (that only depends on E and E') by: 

(29) E" := diag( 7 i/ai,..., 'Jd/a'd ), with := — V 1, k = 1,..., d. 

Pk 

Remark that, by definition of E", one has n<E" = n<E' and n>E" = 'Yqtt>E. Further, 
E" — 70 E is strictly positive on 'H < . 

We can now state our main result. 


Theorem 5. Assume p G D(E,u) with E = diag(l/ai,..., 1/ad) with a.i G (0,1 ),i = 
1,..., d, and a\ € (0,1/2) if d = 1. Let E' = diag(l/a' 1 ,..., l/a’ d ), with a[ > 0, i = 1,..., d, 
and 70 defined as in (1281) . If q^yE) < 2, then 


/ gf(t)-E(5f(t)) \ 

l n'yo+q(E')~q(E")/2 J 


te[o,i] d 


^ (^(^))te[o,i] d ! 


as n —> oo, in the Skorohod space H([0, l] rf ), where (W(i)) tg ]Rd is a zero-mean Gaussian 
process with covariances given by 


Cov(W(t),W(s)) = a 2 x | Yl Co v(B 1/2 (t k ),B 1/2 (s k )) 

v k£l< 


n 

. fce/> 


t'k^k 

~2V 


'n- 


iiog^wr 2 n 


k€l = 


2vr| y k ^ 


with Bi /2 a standard Brownian motion on M, is the characteristic function of u, and is 
given in CD- 


This theorem reveals that taking different summing rectangles may lead to different limits, 
under different normalizations. To the best of our knowledge, such a phenomenon has not 
been noticed in the literature until very recently [33i . [34j for a different model. We elaborate 
more this phenomenon of scaling transition in Section [5j 

a 2 _ 

Remark 4. Observe that one can write Cov(lF(t), W(s)) = s) with 


(30) C(t, s) := | FT t k s k \ / |log if(ir>y)\ 
v fc G /> J jRd 

because of the identity [13, Proposition 7.2.8]: 


_ 2 ( yj ( e itkVk - 1) [e iSkVk - 1) 


n 

,fce/< 


dy , 


( 31 ) 


/ 

Jr 


(fity _ — lj 

- 2 7r |^|i +2 jj- - dy = Ch Cov ( B H {t ), B H {s)) , t, s G M, H G (0,1) 
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with 

C H = --- 

HT(2H)am(Hn) 

The rest of the section is devoted to the proof of Theorem [5j Using Proposition [3l we get 

(32) Sjf(t) - E(Sf (t)) = ]T b n,jWh 

jez d 

with b n {t) = (b n j(t)) jeZd € £ 2 (Z d ) and 

(33) b n j(t ) = Qk—j- 

k£[0,n E 't—1 ] 

Recall that (Xj)j £Z d are stationary martingale differences. 

The proof of Theorem [5] is now divided into three steps. The key step is to compute the 
covariance, which is done in Section ECU Then, we proceed with the standard argument to 
show the weak convergence by first establishing finite-dimensional convergence in Section 14.21 
and then the tightness in Section H. 31 The matrices E and E', and thus 70 and E", are fixed 
as in the assumptions of the theorem. 


Remark 5. As we will see below in the proof, essentially we establish an invariance principle 
for linear random field {Xj)j eZ d with 

Xj = J2qiXj,j€ z d , 


where (X?) ieZ d are stationary martingale-difference innovations, and (qi)i G z d are real Fourier 
coefficients of certain function Q(t). This is a standard framework to obtain linear random 
fields in the literature, and we comment briefly on connections between our results and others, 
(i) First, the same invariance principle should hold if the innovations are replaced by other 


weakly dependent random fields (weakly dependent in the sense of e.g. i, m 0). 


m 


(ii) 


(iii) 


These results can be viewed as generalizations of the seminal work of Davydov 
on invariance principles for linear processes. 

Second, from the modeling point of view, the specific choices of Q(t ) (in terms of 
fi € D(E,u)) and hence (qj)j & z d are new. However, although our assumption on 
Q(t ) is very general, not all possible operator-scaling Gaussian random fields can 
show up in the limit; in particular the Hammond-Sheffield model in high dimensions 
does not scale to fractional Brownian sheets except for a few cases in terms of Hurst 
indices shown in Proposition [ 6 j The aforementioned results [H, 45;] all include 

linear random-field models scaling to fractional Brownian sheets, for flexible choices 
of Hurst indices. 

At last, when the innovation random fields exhibit strong dependence, the limiting 
object could be more complicated ( 0 )- 


4.1. Covariances. 

From (fB2T) . we obtain for f,s£ [0, l] rf , 

Co v(S%\t),S%'(s)) = <Jx(b n (t), b n (s)), 

where, ( b n (t ), b n (s )) := b n ^(t)b nt k(s). The asymptotic behavior of the covariances are 

given in the following lemma where u n ~ v n stands for limn^oo u n /v n = 1 . 

71—>■ OO 
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Lemma 9. For all t,s S [0, l] d , 

a 2 x (b n (t),b n (s )) ~ n 2 ^ +2 ^-^Cov(W(t),W(s)). 


n—>oo 

Proof. Define for m € N and i£l, 


(34) 


A n{x) = 


= > e ilx = 


1=0 


e i(m+l)x __ | 

e ix - 1 ’ 


and for x € M d , the trigonometric polynomial 


K n (t,x ) = Y t je = 


fc=i 


where |_-J stands for the integer part. Recall that since 

Q(x) = £ 


the sequence b n (t) (defined in (l33l) l is obtained by the convolution product of the Fourier 
coefficients of K n (t, •) and Q with Q(x) = Q(—x) since {qj)j & z d * s a rea l sequence. It follows 
that b n k(t) is the fc-th Fourier coefficient of QK n (t, •). Therefore, using Bessel-Parseval 
identity, we get 

{bnif)i bn(s')') — TnZFd I Q(x')K n (t^x')Q(x')K n [s^x')dx 


(35) 

where 


(2tt) 

1 

{2-K) d 

n -FE") 


[ \Q(x)\ 2 l[D 

/ 

J n 


I 1/a' -.iCdOA 1/a' .X x k)dx 


( 2 ^) Jn E " [— 7 T, 7 r] 


<F n (y,i,s)dy, 


*n(v,t,s):= Q(n E "y)\U D l n^U k -ifr 7 fcK 2 /fc)^, n i/a' st _i|(n 


fc=i 


[n ' fcs/c-lj 


According to Lemma [2] and the ^-homogeneity of log?/', one has 


n 


-270 


Q(n E y) =n 270 g(n E y) log ip(n ~ 1oE n *- E ' roE ^y) 


-2 


g(n E y) log y>(n ^ 7 o£ ^y) 


-2 


Thus, 


lim n 270 

n—>■ 00 


Q(n ^'y) = |logV»( 7 r>y)| 


-2 


because E" — 70 FI is null on 'H> and strictly positive on H < and y(0) = 1. Further, for all 
n 6 N*, y € n E " [—tt, n] d , 

2 

(36) 


n 


-270 


Q(n E "y) < max |y(*)| 2 sup | log^(z + 7r>y)| 

cce[- 7 T, 7 r] d z&i< 


-2 
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Now, remark that for all t £ [0,1] and y £ M, 


lim n D [nt _u(n 7 y) = < 

n—>oo L J 


e ity - 1 
iy 


and if \n 7 y| < 7 r, then 


n 


-i 


Di n t-n(n 7 y)| = 


sin(|_ntjn 7 y/ 2 ) 


nsin(n 7 y/ 2 ) 


< 


7 T min 1 , 
7r 

k 2 


if 7 = 1 

if 7 > 1 


\y\ 


if 7 = 1 
if 7 > 1 


where we have used that ^\x\ < | sin(x)| < |x| A 1 for x £ [—7t/2,7t/ 2] and that |t| < 1. Since 
7 fc > 1 if and only if k £ />, we infer 


(37) <S> n (y,t,s) ~ n 270 + 2 9( E ') |logV»( 7 r>y)| 2 [ JJ ] ( JJ —- 


,fce/> 




as n -A oo and for all t, s £ [0, l] d , 
(38) 


n 270 2 g(i?,) |$n(y,i,s)| < max | 5 r(*)| 2 /i(y), 

£cG[— 7r,7rl d 


with 

(39) 


h(y) := sup | log ip(x + 7r>y)| - 2 n min 11 , -—— 

k£l<^ ^ 


Applying the dominated convergence theorem to (1351) , (157)1 and (1381) and using (1301) , to show 
the desired result it remains to prove that h is integrable on K^. 

Formally, write 


h(y)dy = / / h(y) d\ < <g> X>(y) 


n mm {'■ li} iX <M J H> - M P I+v)r 2 n ““I 1 . dX >Jy)- 


kei K 

By Fubini’s theorem, h is integrable over if 
(40) 

and 


J min < 1 , 

fcG/< 


1 


\Vk? 


d\<(y) < oo 


(41) [ h(y)d\>(y)= J sup |log^(® + y)| 2 TT min { 1 , —^ 1 dX>{y) < oo 

Ju> JH> x£'H < l ml ) 


The integrability condition (1401) is obvious. For m , let us remark that the function y £ 
H> t-A inf £Cg -^ < | logip(x + y)\ is ( 7 r>.^-homogeneous and since q{'n>E) > 2, by Lemma [3J the 
function y £ 77 > ga sup 3 , gW< | logip(x + y)\~ 2 is locally integrable on T~L> with respect to A>. 
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Together with the fact that sup^g^ | log ij}{x + y)\ 2 is bounded by 1 for -K = y large enough, 
this shows that 



1 {\\n>y\\ n>E <i}H y ) d ^>(y) < 


with the definition of || • \\ n E given in (HID . Moreover, 



1 {h>y\W >E >i} h (y) d\>(y) 


< 



l{||y||_ B >i} sup \logtp(x + y)\ 2 d\ > (y ) 

n> xG'Hc 



d\=(y). 


The second integral is clearly finite. For the first one, since y € %> H > inf £Cg ^ < | logger + y)\ 
is (7r>£')-homogeneous and q(ir > E) < 2, one has 



lflli/11 B >i} sup I log ip(x + y)\ 2 d\ > {y) 

n> x£'H < 



r q(n>E)-3 



sup I log ^(* + 0)1 2 da n> E{0) < oo, 
a;EW< 


where S 1t> e is the unit sphere of PL > with respect to || • \\ n> E and a 7r> E is the Radon measure on 
Si T> e such that d\ > = r q ^ >E ^~ 1 drda 7T> E- This shows that (1TO holds and thus the function 
h in (HTT) is integrable over M. d . □ 


4.2. Finite-dimensional convergence. 

We start by showing that the coefficients b n j(t) defined in (|33l) satisfy the condition (1201) 
of Theorem [3] in the following lemma. 


Lemma 10. For all t € (0, l] (i and all q = 1, ..., d, 


lim 

n—>■ oo 


1 

1W 


Y - b 2 n}j+eq {t)\ = 0 

j€jj d 


and (E0|) holds. 


Proof. Fix £ € {1,..., d} and t £ (0,1]^ be fixed. Using Cauchy-Schwarz inequality, 


Y 1- b l,j+ez(t)\ < Y (Vi(*) - b nJ+ee (t)) 2 ) 2||6 n (t)||. 

\j £ 7 j d 

So, it is enough to show that 

Y (W*) “ b n,j+e e (t)) 2 = 0(||bn(t)|| 2 )- 


But, we have 


bnj+eiiP) ^ ^ Qk—j ^ ^ Qk—j—ef 


ke[0,n E t-1] 
with k^= \n 1 / a £t£\ — 1 


k£[0,n E t— 1] 
with k£=0 
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Thus, 


/ 


(b n j(t) 


D n,j+e e 


( t )) 2 < 2 £ 


\ 


Let £ > 0. Using Lemma El we get 

/ \ 2 


lin ^ p lMhP ? 


ke[0,n E ' t—l] 

\ with ki=0 ) 


< limsup 

n^-oo ||6 n (t)|p . 


= lim sup 

n—> oo 


( ik-j 

k£[0,n E t— 1] 

\ with kg =0 / 

/ 

1 S 2 ^ k -3 

fcG[0,n s 't-l] 

\with fcf <£'n}^ ol l t( — 1 ) 

\\b n (ti, ■ ■ ■ ,te-i,£te,te + i, ..., || 2 


\ 


IM*)II 2 

U(U,... ,te-i,£te,tg + i,... ,td) 

~ V(t) 5 

where V(t) := C(t,t ) with the covariance function C(-,-) defined in (1301) . We conclude the 
proof of the lemma using that, for any t E (0, l] rf , 

V(ti,... ,t£-i,ete,te+i,... ,t d ) -)• 0, as e -A 0. 

The fact that (1201) holds is a consequence of Lemma [HJ □ 

To prove the finite-dimensional convergence, we use the Cramer-Wold device. Let m E N, 
[0, l] d , Ai,..., A m el, and consider S { „ m) = YX=i A kS^ {t k ). One has 

s<”> - E(S«) = XI dnjX], 
j&, d 

where d n j := EfcLi ^kb n ,j(tk) and Var (S'n”') = ||d n || 2 Var(Xg). Using Lemma El we get 

m m n 2 70 +2g(E')-g(E") ™ m _ 

11 d n 11 = 2_^ ^ k ^i(bn(tk), b n (te)) - (2-kY -/ . / . ^k^eC(tk,tj), 

k= 1 1=\ n ^°° ' n ' k=i l= l 

where C is defined in (130]) . 

If EfcliZXi h^eC(t k ,te) = 0, then nl0+q(E >]_ q(E » )/2 (S% n) - E(sl m) )) converges to 0 in 
L 2 . If ELi ET= 1 *kW(t k ,t e ) > 0, we get that for each k = 1,..., m, 

IIMDII 2 ~ KII 2 tv 

n ^°° Efc=! £*=i A fc A^C(tfc, t £ ) 


Thus, since the b n j(tk) satisfy (1201) . 


sup |d n ,j| < y^Afc sup |b» ; j(tfc)| = ^ A fc o(||6 n (tfc)||) = o(||d„ 

J fc=i J 


fc=i 


This proves that (1201) also holds for the d n j and Theorem [3] applies to . We thus proved 
the finite-dimensional convergence. 
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4.3. Tightness. 

To prove the tightness, by Bickel and Wichura [2'], following [III] and [22], it is enough to 
show that for some p > 0 there exist 7 > 1 and C > 0 such that for all t = (f 1 ,..., id) £ [0, l] d , 


E 


Sf(t)-E(Sf(t)) 


re 


7o+g(£')-^ 


< 


cm 

h-i 


Recall from the equation (|35l) that for all t £ [0, l] rf , we have 
n ~i{ E ") 


\\bn{t)\\ Z = 


to \d [ Q( n E "y) fnk ii( n lk/a ' k yk) 1 


dr/. 


For any 5 £ (0,1), observe that |sin 2 (x)| = | sin 1 ~ 5 (x)|| sin 1 + 5 (x)| < min{|x| 1_ ‘ 5 , |x| 2 } for all 
x, and | sin(x)| > = |x| for x £ [— 7 t/ 2 , 7 t/ 2 ]. Then, for all n and y such that \ny\ < ir and all 
t £ [ 0 , 1 ], one has 


n 2 |Z? Lnt _ 1J (n i y )| 2 = 


— 1 „.\ |2 


sin^Lnfj^) 
n 2 sin 2 (£) 

< min 


.2 j . 1—5 _2 


7T t 


IT 


and thus, 


n 


-2 


£>1^-11 ( n 7 ^)l 2 ^ < 


2 1_<5 |y| 1+5 ’ 4 


7T 1 r . 
ri 1_d min 


< 


7T 


2 1 - 5 


i 1 S min 


\y\ 


1 

1+5’ 1 


2 l-6 
7 r 2 


iz/i 


l+<5 


, 1 if 7 = 1 
if 7 > 1 


Recalling that 7 fc/a], > 1 if and only if A: £ I > , together with (1361) . this shows that there 
exists a constant C > 0 such that 


n - 2 10 - 2 q ( E ')+ q ( E")^ bn ^2 



One can show that this last integral is finite by proceeding exactly as we did to show the 
integrability of the function h in (1391) . The important point is that 1 + <5 > 1 to guarantee 
the integrability of j^pR a t infinity. Hence, for a new constant C' > 0, 




sc'II'l 

3= 1 


—<5 


n -2T0-2 9 (B')+<7(S")||b n (t)||2 < C' 
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Let p > 2. Using Burkholder inequality and the preceding inequality, there exists a constant 


c p > 0 such that 


E 


5f(t)-E(5f(t)) 


n 70 +q(E')-q(E")/2 


— C P^ ( 


b ls lt > 


xf 


n 


2 70 +2 q(E')-q(E")^3 


— Cp 


l|bn(t)|| S 


n 2 7o +2 q(E')-q(E") 
d, 

{l-S)p/2 


< c p c ,p / 2 n tj 

j =i 

which gives the tightness by choosing 5 > 1 — -. 

5. Properties of the limit field 

In this section we focus on the zero-mean Gaussian random field (W(t)) tgK d arising in the 
limit in Theorem [5j Recall that this random field depends on both E and E'. 

5.1. Increments. 

We may consider a harmonizable representation of W, defined on the whole space by 


W(t) = u Y 



|logV’( 7 T>y)l M(dy),vte 


with ax given in (1211) . and At is a centered complex-valued Gaussian measure on with 
Lebesgue control measure (see 0). The harmonizable representation shows that the random 
field has stationary rectangular increments. In the sequel we let (ei,..., e^) denote the 
canonical basis of M rf . Rectangular increments of W are defined for s < t by 

^([M]) = ^2 (—l) d+|e|l IU(si+ £l(tl -Sl),...,S d + £d(td-Sd)) 

ee{o,i} d 

where |e|i = £± + ...+ Ed and corresponds to the directional increment of step 5 € M in 
direction j for 1 < j < d, defined by 

A { pW(t) = W(t + Sej) - W{t). 

A direct consequence of Theorem [5] are the following properties of the random field W. 

Proposition 4. The random field W satisfies the following properties: 

(i) stationary rectangular increments: for any fixed s € 

(W([s, t])) s<t f = (W([ 0, t - *])).<* = (W(t - s)) s<t ; 

(ii) (E', H)-operator-scaling property: for all A > 0 

{W{\ E 't)) tmd f = (X H W(t)) tmd , 

_ q(E») 

2 • 


with H = 70 + q{E') 









INVARIANCE PRINCIPLES FOR OPERATOR-SCALING GAUSSIAN RANDOM FIELDS 


27 


We can say more about the directional increments A ^W(t). First of all, as a special case 
of Proposition U @), W(t ) viewed as a process indexed by tj £ M has stationary increments. 
Moreover, simple dependence properties in the directions corresponding to /> and /<, if 
not empty, are given below. Following ideas from [33;, Definition 2.2] we state the following 
proposition. Recall that |/>| < 1. 

Proposition 5. The random field W satisfies the following properties: 

(i) When /> = {j}, the random field W has invariant increments in the direction ey. for 

all h, 6 £ M, t £ M. d , we have A^W(t + hefi = A ^pW{t). 

(ii) When /< 0, the random field W has independent increments in any direction e,j with 

j £ /<: for all 5 > 0,i £ M d , A^pW(t) is independent from W{t). 

Proof. Let (e.j) L denote the subspace of M. d orthogonal to e r Let T( ej )J- and A( c ^_l denote 
the corresponding projection and Lebesgue measure, respectively. First, let us simply remark 
that for /> = {j}, 5 £ M, and t £ M d , 

(42) A®W{t)=6W('K {e . ) ±(t) + e j ), 

which does not depend on tj. The desired statement then follows. For the second statement, 
when j £ /<, 


(43) Afw(t) = a x [ []4 

<kEl> 


/ 

J R- 


ytjyj fySyj — ^) 


Wj 


n 

K kel<;k^j 


Wk 


|log^(7r>y)| l M(dy). 


Therefore 


Cov(A fw(t),W(t)) = C e 


,(t) [ 

Jw 


( e i8y i - 1) (1 - e*bw) 


\Vj? 


dyj, 


with 


c ej (t)=(j 2 x [ n 4 

\kei> 
,(i) 


L n 


e itkUk _ \ 


Wk 


\\ogfi{-K>y)\ d\ {ej) ±(y). 


Hence, Cov(Ay j W(t), W(t)) = 2irC ej (t)Cov(B 1 / 2 (t j +6) - B 1 / 2 (t j ), B 1/2 (tj)), with B l/2 
a standard Brownian motion on M. By independent increments of Bi/ 2 , we obtain that 

Cov(A^LF(t), W(t)) = 0 for <5 > 0. Since W is a Gaussian field, we conclude that ApW{t) 
is independent from W(t). □ 

Let us quote that our definitions of invariant and independent increments are not the ones 
used in [33, Definition 2.2], However we remark that invariant increments in the direction ej 
lead to invariant rectangular increments in the sense that, for all 5 £ R, and s < t 


W([s + 5ej,t + 5ej])=W([s,t}). 


This follows from the fact that 


W([a + Seyt + Sej]) = A« 8l A^_ S2 ... A%L ad W(a + Sefi. 


(d) 
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Indeed, computing first A [^_ s .W(s + dej) = Awe obtain the desired result. 

When the increments are either invariant or independent in at least one direction, we say 
that W has degenerate increments. Otherwise, we say that W has non-degenerate increments. 

Example 4. When d = 2, choosing E' = diag(l,/3) for > 0 as in [33] we obtain that 
|/ = | = 2 if and only if p\ = p 2 , that is /3 = It follows that for j3 one has |/=| = 1 

and W has either independent or invariant increments in the orthogonal direction. However, 
when /3 = ^ we get 

j- / it^yk _ ] \ 

W(t) = a x / IT -—:- |logV’(y)r 1 A4(dy),Vt G M 2 . 

In this case, W has non-degenerate increments. Recall that all possible non-critical cases in 
d = 2 have been provided in Theorem [2] in introduction. 

More generally for d > 2 we can state the following scaling-transition property. 

Corollary 1. The random field (Xj)j £Z d, defined in Section Wf^. exhibits a scaling-transition 
in the sense that 

(i) If there exists c > 0 such that E' = cE, then W has non-degenerate increments; 

(ii) Otherwise, W has degenerate increments. That is, there exists at least one direction in 
which the increments of the limit random field are either invariant or independent. 


In the sequel, we need to control the variance of the directional increments. By Proposi¬ 
tion [5] for all u G M d , 5 G M, 

(44) Var(A .£V(t*)) = 5 2 Var(W(7r (e , } x (u) + ej )),j G />, 
and 

(45) Var(A fw{u)) = |5| Var(W( 7 r (e . ) x (u) + ej)),j G /<. 

The control for j G 1= is a little more involved, as summarized in the following lemma. 

Lemma 11. There exist some constants C such that for all u G [—1, l] rf ,<5 G 1, j G I=, the 
following inequalities hold. 

(a) If |/>| = 1 or /> = 0 and aj < 1/2, 

(46) Var(A^'*IT(u)) < C\5\ 2 ^ j with (5j = aj ^1 — _|_ I. 

(b) If /> = 0,<x,- = 1/2, then 


(47) 


Var(A^V(u)) < Cmax(^, |<5| 2 ^) for all Hj G (0,1). 


2 | e|2 H, 


(c) If /> = 0, aj > 1/2, then 

(48) Var(A^V(u)) < C5 2 . 

Proof. Recall fpTlT) . For j G /=, for all u G [—1, l] d and 6 G M, 

Va I (A«Vw) = L n «.| j e “ u ■ 1,2 


fce/> 


l Vj 


fj(Vj)dyj , 
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with 


=J n 


3 ' k£l<\k^j 


e ^kUk _ i 


Wk 


I log V’( 7r >( 2 / + U j e j ))\ d\( ej) ±{y) 


This is a locally integrable function over M for all values of ay € (0,1) due to the fact that 
| log V , ( 7r >2/)I is a vr>£ , -homogeneous function, q(ir>E) > 2, and Lemma [3l Furthermore, by 
^-homogeneity and polar coordinate x = t(x) e 0(x), 

llog^aOl- 1 = 1 ' (I log^(7T>a?)l + \xjn~ 1 


T { x )\ log '4’( 7r > 6(x))\ + T(x)\dj(x)\ a j 
r( *)|log^(0(®))| 

< ci(|logV'(7r>cc)| + 

with ci = max 0eSB (|logV>(vr>6»)| + foyfo-fo/l \ogi/)(0)\. Thus, 

2 


(|log^(7r>*)| + \xjp) 


OL A \ — 1 


fj(vj) < c i [ 


J ' k£l<\k^j 


glUkUk _ 2 


Wk 


(I log V>(tt> 2 /)I + I yj\ aj ) 2 d\ {ej) ±(y) 


= ci 


n 

. fce/<;Mi 


gLUkVk _ ^ 


Wk 


dyk jJ H (l 1 °g'0( 7r >3/)l + \yj\ aj ) 2 d\>(y), 


and the second integral in the right-hand side above equals 

, -2 


'Hi 


foil 2aj {\ lo g^((fojT J ) E n>y )| + l) d\>(y) 


= \ yj \f (| log ip( 7r > y)\ + l)- 2 d\>(y) =: foi 
JH> 


C2 


with /3j = ctj( 1 — q(n > E)/2) + 1/2. We have thus obtained 

fj(yj) < c 3 fo j |-^ +1 with c 3 = c 2 c 2 (2vra fc ). 

fce/<;fc^i 

Recall that |/>| < 1. 

(faj) In case that |/>| = 1, q(ir > E) > 1 and thus ,8j < 1. Therefore by the above calculation 
and (13T1) . 


(49) Var(A^W(u^)) < cr\c 3 


e iS Vj _ 1 


Wj 


foil 2 ^ j+1 dyj = a 2 x C3C Pj \5\ 2f> E 


In case that |/>| = 0, / 3j = ay + 1/2. If ay < 1/2, then the same bound (|49l) holds. 
(Ep If ay = 1/2, then for any Hj £ (0,1), 


r e iSyj _ x 2 r 

/ ——— fj(yj)dyj <t 2 fj(yj)dyj + c 3 

Jr Wj J\yj\<X 


2Vj\>^ 


e i5 Vj _ 1 


z 2/i 


foi|- 2 ^ +1 ^i 


< C4max(<5 2 , |(5| 2Hj ), 
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with 


ttef—1,1] i J/ e .u , , 


C 4 = max 


e iUkVk _ i 


iu k 


|log^( 7 T>j/)| 2 d\ {ej) ±(y) + c 3 C Hj - 


Therefore, 


Var ^ W(u ^< 0/^04 max(<5 2 , \8\ 2Hj ). 

(cf) At last, if ay > 1 / 2 , then (3j > 1 , the function fj is integrable on M and 

e iS Vj _ 11 2 


/ 

JR 


Wj 


fj(Vj)dyj < 5 / fj(yj)dyj. 


It then follows that 

with 


Var (A^'V(w^)) < c 5 d 2 , 


c 5 = a 2 x sup / 


e iUkVk _ \ 


Wk 


|log^(vr>y)| 2 dy. 


□ 


5.2. Fractional Brownian sheets. Here we give a complete characterization of when W 
is a fractional Brownian sheet. Recall that a zero-mean Gaussian random field (V(i)) tgR d is 
a standard fractional Brownian sheet with Hurst index (Hi, ..., Hj) £ (0, l] d if 

1 d 

(50) Cov(X(t),X(s)) = ^ ( \U\ 2Hi + \si\ lHi - |U - Sj| 2 ^) . 

i= 1 

Remark that we include the degenerate case that Hurst index equals 1. 

For the limit random field W, the covariance function can be factorized according to dif¬ 
ferent directions as 

2 

Cov(IF(£),lF(s)) = —Cov(B 1/2 (t fe ),B 1/2 (sfe)) • JJ t k s k -^(t,s), 

' n ' kei< kei> 

with \H(i, s ) only depending on {£&, Sk}kel =, given by 

*(«,»)- /1 log !/■{«)r 2 n — ~ 2 1 ? (e ;r* ~ ~ dx >(v'i■ 

Jn> R. 2»l»l 2 

Recall Ch in (13T1) . 

Proposition 6 . Then W is a fractional Brownian sheet, if and only if \I=\ = 1 .In this case, 
T(t, s) has the following expressions: in case 1= = {j},I > = 0, 

(51) V(t,s) = \ \og^(e j )\~ 2 C Q . + i/ 2 Cov(B aj+ i/ 2 (t j ),B Q . + i /2 (s j ))] 
in case I = = {j}, /> = {k}, 

(52) V(t,s)= f | log ip(y + ej)\~ 2 d\ > (y)CHj Cov(B Hj (tj), B Hj (sj)), 
with Hj = atj( 1 — 1/(2 a*,)) + 1/2. 
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Proof. We first prove the ‘if part’. Suppose 1= = {j}. In the case /> = 0, 


^(M) = 


0 (e*b'%' _ i )( e ^jyj _ i) 

|log^(^ e ,)|- 2 -^-- % 


_ life isjVj _ 1) 

llogiKd %r)S)r 2 -- 5 -r-p- 

Z 7 r li/i I 


2 (e itjyj — l)(e is i y i — 1) 

I log V ; ( e j)l "- o„L...I 2 + 2 Q , - d yp 


2it\yj\ 2+2a i 


Thus by (I3IT) . in case I = = {/},/> = 0, (15T1) follows. In the case I> / 0, 

( e *b%' _ i)( e ^w _ i) 

i iu s wyy v vj^j)' 
m> 


T(t,s) = / / |log^(i/ + i/jCj)|' 

JR 


27 r bil 2 


-d\ > (y)dy j 


Sfep- dA > (!,)c, w 


/ [ | W |- 2 “'|logV’((ter“>) E (!/ + we j ))|- 2(e '‘''" 1)(e “ ,W x) 

JR JH> 

f of - 1) 

/ llog^ + e,)!- 2 ^^) / 1 j 

J-H> JR 


27r|y J | 2+2a J _a ^( 7r> - E ) 


IH> J M 

That is, in case /= = {j},/> / 0, for Hj = Oy(l — q(ir > E)/2) + 1/2, (|52l) follows. 

Next, we prove the ‘only if part’. Suppose W is a fractional Brownian sheet with Hurst 
indices Hi, ..., H ( j. From Proposition 0] IF is also (E 1 , i/)-operator-scaling with H = 70 + 
q(E') — q{E")/2. Then, it follows that 


— + • • • H-- — 70 + q(E') — q(E")/ 2, 


a 


a„ 


or equivalently 
(53) 


fce/< k kei> k \ keiy K J 


We consider the variance. By the assumption that IF is a fractional Brownian sheet, and the 
fact that IF has stationary directional increments, for all j € {1,..., d}, for all 5 E M, 

(54) Var(A^IFH) = |<5| 2 ^ Var(IF( 7 r (e . > x(w) + e,)). 

Recall that |/>| < 1. We first consider the case /> = 0. In this case, 

• for k £ /<, comparing ((Ml) and ([751) yields H *. = 1 / 2 , 

• for k € I=,ak < 1/2, comparing (1571) and (061) yields Hf = otk + 1 / 2 , 

• for /c € I=,&k = I/ 2 , comparing (1571) and (071) yields Hf. = 1 , 

• for A; € /=,«& > 1 / 2 , comparing (1571) and (051) yields Hk = 1. 

Then (l53l) becomes 

E E 7. = 7C 

k£l=,a k > 1/2 k£l=,a k < 1/2 

Since a*. < 1 , it then follows that |/=| = 1. Similarly, in the case /> 7 ^ 0, say /> = {1}, it 
follows from comparing the corresponding inequalities that 

• Hi = 1 , 

• for k £ /<, Hk = 1 / 2 , 
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• for k £ I=, H k = a k ( 1 - l/{2a 1 )) + 1/2. 

Then, (|53l) becomes 

which implies |I = | = 1. □ 

Remark 6. When the limit is a fractional Brownian sheets, in directions corresponding to 
I>, /< (if not empty) and /=, the Hurst indices equals 1, 1/2 and some value in (1/2,1), 
respectively. Thus, W exhibits long-range dependence in the directions corresponding to />. 

As a concrete example, we prove Theorem [2] 

Proof of Theorem [H Case {*)): when a 2 > a 2 ,a 2 £ (0,1/2). In this case, 70 = P 2 = ct 2 /a 2 , 
E" = diag(l/ai, l/a' 2 ), /< = {1},/= = {2}, /3 = a 2 /a' 2 + ~(^- + ^) and = 1/2 are 

straight-forward. Then, by (J5TJ) , H 2 = \ + a 2 and cr 2 = Ch 2 \ log -0(0,1)|~ 2 . 

Case ([HJ/: when > 02,0:2 £ (1/2,1). In this case, 70 = p\ = 1, E" = £’,/> = {2},/ = = {1}, 

/3 = l+ 2^- + ^r-2^ and H 2 = 1 are straight-forward. Then, by ([52]), Idi = i + aq(l - ^) 
and a 2 = C Hl f R | logif(l,y)j~ 2 dy. 

The calculation of cases Iml) and Ini) are similar and thus omitted. One obtains a 2 = 

Ch J logi/’(l,0)|~ 2 for case Im1) and cr 2 = C# 2 | log?/(y, l)|~ 2 dy for case IS)) . □ 

5.3. Sample-path properties. 

We conclude this section by the following general sample-path properties for the random 
field W that is a consequence of Proposition 5.3]. 

Proposition 7. There exists a modification W* ofW on [0,1]^ such that for all e > 0, almost 
surely there exists a finite random variable Z such that for all s,t £ [0, l] d , 

I W*(t) - W»| < Zp(t, s)log(l + p(s, t) _1 ) 1/2+e , 

with 

P(M) = Y, I tj - s il + X] to “ s j ! 1/2 + to “ 

jei> jei< jel= 

where for j £ /=, 

(a) Hj = otj (1 — q('K > E)/2) + 1/2 if either |/>| = 1 or /> = 0 and a 3 < 1/2, 

(b) Hj can take any value in ( 0,1) if I> = 0 and a.j = 1/2, and 

(c) Hj = 1 if I > = 0 and Oj > 1/2. 

Proof. Let us consider the diagonal matrix with entries corresponding to 1 for j £ />, 2 
for j £ and 1 /Hj for j £ /=. Let te»> be the radial part with respect to E'" according 
to |5, Equation (9)]. Let us quote that since t 1 -+ p(0,t ) is E'" homogeneous and strictly 
positive on \ {0}, following ideas of Clausel and Vedel Theorem 3.2], the function 
t I-+ p(0,t)/T E III (£) is continuous and strictly positive on the compact set Sf" 1 ■ It follows that 
we may find C, C' > 0 such that for all t £ M d , 

Ct E '"(£) + p(0,t) < C't E "i( t). 
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Therefore, by Proposition 5.3] (with (5 = 0), to show Proposition[7]we prove for t, s £ [0, l} d 
that 

(55) y/E((W(t) - fP(s)) 2 ) = ^Var (W(t)-W(s)) < Cp(s, t). 

For t, s € [0, l] rf , considering as in 0, the sequence (rt^^)o<j<d defined by u = s and 
w (i+i) = U (J) + (tj — Sj)ej for 0 < j < d — 1, we get W(t ) — W{s) = A^_ g ^W(u^). 

Hence 

d 

(W(t)-W(s))<J2 

3 =1 

Now to obtain (1551) . it suffices to apply the bounds on the directional increments established 
in Lemma EH Observe that in the case j £/=,/> = 0, since 6 = tj — Sj € [—1,1], the right- 
hand side of ff7l) becomes C\5\ 2H i. The details are omitted. The proof is thus completed. □ 

Let us quote that we probably could improve this result. Actually, following [Ifij], it is 
sufficient to get a similar lower bound on the variance on [e, l] d to establish condition (Ci), 
from which Theorem 4.2 follows, saying that the inequality is true for e = 0 and Z has finite 
moments of any order. 
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